ON THE ASYMPTOTIC DISTRIBUTION OF ZEROS 
OF MODULAR FORMS 

O 

^ ■ ZEEY REDXICE 

C . 

^ , Abstract. We study the distribution of zeros of holomorphic 

modular forms. Assuming the Generalized Riemann Hypothesis 
we show that the zeros of Hecke eigenforms for the modular group 
become equidistributed with respect to the hyperbolic measure on 
the modular domain as the weight grows. 

H 

1. Introduction 

EL Our purpose in this note is to study the limiting distribution of 
zeros of modular forms. We review some definitions: A modular form of 
weight k for SL 2 {'L) is a holomorphic function on the upper half-plane 



(N 



M, transforming as /(gg) = (cz + d) k f{z), for all ( a c b d ) G SL 2 (Z) 
(this forces k to be even), and "holomorphic at the cusp" (see § I2.1JI . 
A form is cuspidal if it vanishes at the cusp. For a modular form of 
^ ■ weight k, let v(f) be the number of inequivalent zeros of / in H, with 

q ■ the convention that a zero at z is counted with weight w(z) inverse to 



the number of elements of iST 2 (Z)/{±_/} fixing z. Then u(f) < k/12. 

For a sequence of modular forms, where we assume that the number 
of inequivalent zeros tends to infinity, we would like to examine the 
manner in which the resulting configuration of zeros is distributed in 
the modular domain SL2(Z)\EI. 

X 

1.2. Example: Eisenstein series. The Eisenstein series are non- 
cuspidal modular forms of weight k > 2, given by the sum 

F.K.C. Rankin and Swinnerton-Dyer showed that all zeros (in the 
fundamental domain) of lie on boundary arc {\z\ = 1}. Moreover, as 
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k — ► oo, the zeros become uniformly distributed on the unit arc. R.A. 
Rankin ^3] gave a similar result for the (cuspidal) Poincare series 

with m > 1 fixed, k^> m. See |T] for yet another such example. 

1.3. Example: higher Weierstrass points. An example in a slightly 
different setting are the higher Weierstrass points on a (compact) Rie- 
mann surface C. Let Sk(C) be the space of holomorphic forms of weight 
k on C, d := dxmSk(C), and let /i, ...,/<; be a basis for Sk(C). The 
Wronskian is defined as 

Wrfctf/,-};*) := det(if 

It is a modular form of weight d(k + d — 1). 

A k-th order Weierstrass point is a zero of Wr^. Equivalently, these 
are points of C where there is a nonzero form which vanishes to order 
> dim S k (C). 

It was conjectured by Bers, and proved by Olsen [TU] that as k — > oo, 
these become dense in C. Mumford (see page 11 of [7]) and Neeman jH] 
showed that in fact these become equidistributed with respect to the 
Arakelov (or Bergmann) measure on C, which arises from the metric 
on C gotten by embedding the curve in its Jacobian and pulling back 
the flat metric. 

1.4. In contrast to these examples, we consider the case of cuspidal 
Hecke eigenforms. Recall that the Hecke operators 

ad=n b mod d 

act on on the space of cusp forms of weight k, commute with each other 
and are self-adjoint with respect to the Petersson inner product 

k dxdy 



(f,g):= / f(z)g(z)y k , 

JSL 2 (l,)\M V 

Thus the space of cusp forms admits a basis consisting of joint eigen- 
functions of all Hecke operators (the Eisenstein series Ej. is also an 
eigenfunction). The cuspidal Hecke eigenfunctions have a simple zero 
at the cusp. For these Hecke eigenforms, we have: 

Theorem 1. Assume the Generalized Riemann Hypothesis (GRH). Let 
{fk} be a sequence of cuspidal Hecke eigenforms, then as k — ► oo their 
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zeros are equidistributed in SX 2 (Z)\HI with respect to the normalized 
hyperbolic measure dV(z) = |— p-. 

Equidistribution here means that for any nice compact subset Q C 
SX 2 (Z) \H, the proportion of (inequivalent) zeros of such a Hecke eigen- 
form which lie in Q is asymptotically the relative area of Q: 



1.5. Theorem n is the outgrowth of some ideas from the theory of 
"quantum chaos" . The mechanism is that equidistribution of the zeros 
of a sequence of forms fk is implied by equidistribution of the "masses" 
y k \fk{z)\ 2 dV(z) of the forms. This idea was discovered by Nonnen- 
macher and Voros in the context of quantum maps. Around the 
same time a very general result of this sort was obtained by Shiffman 
and Zelditch ^H] for zeros of sections of high powers of a positive holo- 
morphic hermitian line bundle over any compact complex manifold ([9 
deal with curves of genus one). All these are in a compact setting. The 
analogous result in our case is 

Theorem 2. Suppose that {f^} is a sequence of L 2 -normalized cusp 
forms (fk of weight k) for which the bulk of zeros lie in the fundamental 
domain: v{fk) ~ k/12. Assume that for some c > 0, we have 



Then the zeros of fk are equidistributed with respect to dV(z). 

Here — > denotes weak convergence when we test against compactly 
supported functions. 

The proof of Theorem |21 given in Section |SJ follows closely the argu- 
ment of [HI EE], with care taken to handle the complications due to the 
lack of compactness of the modular domain. We require the hypothe- 
sis which is weaker than equidistribution of the masses as it allows 
"leakage" of some of the mass at the cusp (which cannot happen in a 
compact setting). This slightly weaker version of equidistribution of 
masses is what Lindenstrauss [5] proved unconditionally for the anal- 
ogous case of Maass forms, though he cannot exclude c = 0. In the 
holomorphic case we do not have an unconditional proof of (Q). How- 
ever it has been known for some time to follow (with c = 1) from GRH 
[TT)1 12T] . This is sketched in section HJ 
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y k \f k (z)\ 2 dV(z) 
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2. Potential theory on r\H 

2.1. Preliminaries. We review some definitions: A modular form of 
weight k for SX 2 (Z) is a holomorphic function on the upper half-plane 
H, transforming as 

/(^) = (cz + <*)*/(*), for all (l G SL 2 (Z) 

(this forces k to be even), and holomorphic at the cusp. This means the 
following: Since f(z) is periodic, it can be expressed as a holomorphic 
function f(q) of q — e 2mz in the punctured disk < \q\ < 1. The 
requirement to be holomorphic at the cusp means that it extends to a 
holomorphic function at q — 0. A form is cuspidal if it vanishes at the 
cusp, that is /(0) = 0. The order of vanishing at the cusp ordoo(/) is 
defined as the order of vanishing at q = of f{q). 

We denote by V = SX 2 (Z)/{±/} and speak interchangeably about 
modular forms for T and for SX 2 (Z). For each z G H let T 2 = {7 G 7 : 
72 = 2} be the stabilizer in T of z, and set 

zeYi 

j g rV 

otherwise 




:= 77^ = < 1/3, z G Te 2 ^ 3 

#1 2 



For a modular form of weight k, let 

z£SL 2 (Z)\H 

be the (weighted) number of T-inequivalent zeros of / in HI. Then 
< t|, in fact + ord oc (/) = A Note that = only for 
powers of the modular discriminant. 

2.2. For a smooth compactly supported function <fi G C£°(M) on the 
upper half-plane H set 

f^z) :=^0( 7 ^)Gcr(r\e) 

Note that 

/ F ( p(z) ( ^^- = I (f)(z)dxdy 
Jr\m U Jm 

Let dV(z) = vol (r\ H ) be the normalized hyperbolic measure on 
the quotient T\H. 
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Lemma 2.1. Let f be a (weakly holomorphic) modular form of weight 
k for T and let {zj} be a set of T-inequivalent zeros of f in HI. Then 



(2) ^ w (^.)F^-) = *^ffi/ F,(z)dV( 
i 47r Jt\w 



J- [ \og(y k / 2 \f(z)\)Acf>(z)dxdy 



+ 2, 

Proof. Let S = {z e EI : /(z) = 0} be the set of zeros of the form /. 
It is a discrete set of points, which is stable under V (^S = S for all 
7 G T). For any such S we have 

(3) £ u,^^ = 
Indeed, 

E^)= E I» 
= E E *w 

2er\5 7 er/r z 

= E #^E^> 
= E 

as required. 

We recall that the electrostatic potential for a point charge in the 
plane is ^ log |z|, that is for e C£°(C) we have 

^log|z|A0(z)^y = 0(O) 

Consequently, for the holomorphic function f(z) on H we have 

(4) / ±]ag\f(z)\&Kz)dxdy= £ (f>(s) 

for all <fi E C£°(H) (with multiple zeros repeated). Thus 



= 1- J log\ f(z)\A(i>(z)dxdy 



.l-Jjogy k / 2 A ( f>(z)dxdy + ^- \ \og{y k / 2 \f{z)\)A(t>{z)dxdy 
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The first term above is transformed via integration by parts into 

-i- / A(logy k / 2 )<t>(z)dxdy = A / ftz)^ 

2vr J B 4vr J m y 2 

Note that 

/ <Kz) AA = / F,(z)^ = vol(r\H) / F,(z)dV(z) 
so that we get 

~ [ A(\o gy k / 2 )4>(z)dxd y = fc VQl( ^ W) I F^z)dV(z) 
as required. □ 



Remark: We may reformulate this in T-invariant form as follows: 

a 2 , a 2 

dx 2 dy 2 



Instead of the Euclidean Laplacian A = + we use the hyperbolic 



Laplacian 

d 2 d 2 



£ = y' 



dx 2 dy 2 



which is T-invariant, and instead of \f{z)\ we use y k \f{z)\ 2 which is 
T-invariant. Then 

(5) ^^) = fe!ffl/ F 4l {z)dV{z) 



^ voi(r\H) 



log(// 2 |/(,)|)£F (z)^ 
2tt J r \ m y z 



To derive © from Lemma \'2.1\ we transform the second term in (J2J) by 
noting that since is T - invariant, we have 



A / iog(// 2 |/WI)^(^)AA = 1/ iog(^ 2 |/(,)|)^(/:0)( 7 z)^ 

= A/ log(// 2 |/(,)|)^(.)^ 
27T Jr\ H r 



vol(T\H) 



Mv fc/ V(*)l)*/*(*W*) 

r\H 



2tt 

This proves (JSJ) once we note that F C( j, = CF^ since C is T- invariant. 
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2.3. Background on sub-harmonic functions. Let O C C = f 2 be 

an open connected set. An upper semi-continuous 1 function u : Q — > 
[— oo, oo) is sub-harmonic if u is not identically — oo and for all zGfl, 

1 f°° 

u{z) < — / u{z + re ie )dO 
2vr Jo 

whenever the disc {w : |z — w\ < r} is contained in fi. 

Note: m is sub-harmonic -v^ Au > as a distribution. 

A fundamental example of a sub-harmonic function is log \f(z) | where 
f(z) is holomorphic in Q. Moreover we have 

Alog|/(z)|=27r 5 ( z ~ z i) 

/(*j)=0 

We denote the space of subharmonic functions on Q by S7i{Q). 

A basic compactness property of sub-harmonic functions is (c.f. [2|, 
Thm 4.1.9): 

Lemma 2.2 (Compactness property). If {uj} C S7i(Q) are locally 
uniformly upper bounded (that is for all compact K C Q there is ck so 
that Uj(z) < ck for all z G K and all j) then either 

(1) Uj — > — oo uniformly on compacta, OR 

(2) There is a subsequence Uj k which converges weakly to some u G 

snip). 

(3) Moreover, in this case lim sup Uj < u and lim sup Uj = u almost 
everywhere. 

The following is known as "Hartogs' Lemma" (^H] Thm 3.4.3): 

Lemma 2.3. // {uj} C SH(Q) are locally uniformly bounded above, 
and there is a continuous function <fi : Q —>■ R such that lim sup Uj < 
4>, thenmeLx(uj, <p) — > <fi locally uniformly on Q as j — > oo. 

3. Ergodicity of eigenfunctions implies equidistribution 

OF ZEROS 

3.1. Proof of Theorem [2J Assume that there is some c > so that 
y k \fk\ 2 dV — > cdV for a sequence of k — > oo. We will show that nec- 
essarily the zeros {zj} of fk become equidistributed relative to dV as 
k — >• oo. We need to show that for all F G C£°(r\H) we have 

J—J2 w ( z i) F ( z j)~ [ F(z)dV(z), fc^oo 

V\Jk) j JT\W 

1 u : Q — * [— oo, oo) is upper semi-continuous if for all c £ K, the set {z : u(z) < c} 
is open. 
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It suffices to show this for F of the form F(z) = F^(z) = X] 7 er 4>{l z )i 
4> G C£°(H), since these are dense in C^°(r\EI) with respect to the 
uniform topology. 2 For these we have the fundamental identity of 
Lemma 12.11 



(6) = ~77T / F ^)dV(z)+£ 



£ = 7T^ttt I log(y«\f k (z)f)AcP(z)dxdy 



where 

£ = 

Since we assume that v(fk) ~ k/12, the main term is the desired one 
(the mean value of F^) and we need to show that £ — > as k — > oo, 
i. e. that for all G C c °°(e), 

^log(y fe i/ fe (^)| 2 )A«/-(z)dxdy^O 

or that 

(7) / ilog(|/ fc (^)| 2 )A0(z)^^ / -logyA^dxdy 
Jm & Jm 

Assume that this is false, that is there is some test function <po G 
C^°(EI/r) and a subsequence of the /& (which we will continue to call 
fk for notational ease) for which J7J) fails. We will proceed to derive a 
contradiction. 

We begin by listing some properties of v k ■= \ log These are sub- 
harmonic in H/r. By Proposition El we have y h \fk{z)\ 2 ^ k uniformly 
on compacta and so 

logfc + 0(l) 

(8) t> fc <-logy+ k 

thus 

(9) lim sup v k < — log y 

locally uniformly. Thus the family {v k } C S7i is locally bounded 
above. By the compactness property fLemma 12. 2|) . this gives us two 
possibilities: 

i) Vk — > —oo locally uniformly, OR 

ii) {vk} has a weakly convergent subsequence. 
We will dispose of both possibilities. 



2 To see this it suffices to show that for fixed ball B C H, the space {F^ : 
supp0 C B} is dense in C(n(B)) 1 it : H — >• r\H the projection, which follows by 
applying the Stone- Weierstrass theorem. Note that the uniform closure of the space 
{Fff, : <j> € C£°(H)} is Co(r\H), the continuous functions vanishing at infinity. 
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Option i) On the support of the test function O , we have Vk — > — oo 
uniformly, so that there is some K so that for all k > K and all 
z G supp 0o we have 

v k (z) < -2H 

where H = max{Im2 : z G supp^o}- But then exponentiating we find 

\M 2 < e~ 2kH 

so that for all (f> supported inside supp (fi (and so that supp (p is con- 
tained in a single fundamental domain) 

/ *(*)l/*(*)lV^->o 

Jr\m V 

contradicting the assumption that y k \fk\ 2 dV — ► cdV. 

Option ii) We assume {vp.} has a weakly convergent subsequence, 
which for notational convenience we continue calling {v^}, which then 
converges to some v G SH, and moreover limsupffc < v and these are 
equal almost everywhere. Then by v+logy < almost everywhere. 

Since v k — > v but J m v(z)A(f) (z)dz ^ J„ — logyA(j) (z)dz, we know 
that i) / - log y on a set of positive measure. Since v is upper semi- 
continuous, there is some 5 > so that v < — log y — 5 on some 
nonempty open relatively compact set U. 

Since v = limsupf^ < — logy — 5, by Hartogs' lemma EHfl there is 
K = K(5, U) so that for all k > K, v k < - logy - 5/2 on U. That is 
to say, 

ilog|/ fc | 2 <- logy -6/2 

or 

y k \h\ 2 <e- k& ' 2 

on U, which as in the first option contradicts y k \fk\ 2 dV — > cdV. 

4. Sketch of the proof of Theorem [T] 

As is well understood by now JOIIS], the equidistribution of masses 
(P) (with c = 1!) follows from the Generalized Riemann Hypothesis 
(for certain automorphic L-functions), in fact from a "subconvexity 
estimate" for the central value of certain L-functions. To relate the 
equidistribution of the masses for Hecke eigenforms to GRH, one has 
to examine the behavior of the "periods" 

f g{z)\f k {z)\ 2 y k dV{z) 
JSL 2 (i)\a 

where g is fixed and k — > oo. We need to show that whenever g is 
orthogonal to the constants then the period vanishes as k — ► oo. 
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In the case where g is Eisenstein series the period is the (completed) 
Rankin- Selberg L-function L*(l/2 + it, fk x /*.), the standard convexity 
estimates give a bound of k e (for all e > 0) and vanishing as k — > oo 
follows from GRH (see [13 §4]). 

To handle the case when g is a cuspidal Hecke-Maass form, one 
exploits Watson's formula [21] which relates the period with the central 
value of the triple product L-function: 

L(~,9* fk x f k ) 

To get the decay as k — > oo, it transpires that, again, one has to beat 
the "convexity bound" on the central value, for which one uses GRH. 

In the case of CM-forms 3 , Sarnak ^7] proved that their masses are 
equidistributed. Thus the reasoning above implies (unconditionally) 
that the zeros of CM-forms are equidistributed w.r.t. dV as k — > oo. 

Appendix A. An L°°-bound for cusp forms 

A holomorphic form of weight k for SL2CZ) is bounded, since we 
require that the g-expansion contains no negative powers of q. We 
will need a bound on the supremum with an explicit dependence on 
the weight k. For our purposes any bound on y k ^ 2 \f(z)\ which is sub- 
exponential in k for fixed z suffices. Below we derive a bound of size 
A; 1 / 2 ; since the dimension of the space of cusp forms grows linearly in k, 
this bound is optimal (for arbitrary forms) as far as the fc-dependence. 
The proof is adapted from [3J Lemma A.l] who treat the case of Maass 
forms, though unlike |S a we do not need to control the z-aspect. 

A.l. The incomplete Gamma function. We first recall some prop- 
erties of the incomplete Gamma function T(a,x) defined for x > 
by 

T(a,x) := / e-H a - 



t 

From the definition it is clear that T(a, x) is decreasing in x. For integer 
k we have 

(10) r(fc- l,x) = (k-2)\e~ x J2 

m=0 

What is not completely straightforward is the asymptotic behavior 
of T(a,x) when a,x tend to infinity. We will need the asymptotic: 

(11) r(jfc-l,A;)~-r(Jfe-l), k ^00 



x 
rn! 



3 These only exist for certain congruence subgroups of the modular group 
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By (ITHjl . the relation (JTTJ) is equivalent to 

(12) e -*y* * 

m=0 

This is close to a conjecture of Ramanujan ^TJ EH |2U| that 

1 

~ 2 

for some § < = 0(&) < |; he showed [TJ] that 0(oo) = 1/3. We need 
the weaker asymptotic (|T2*|) or equivalently (by Stirling's formula), that 



(13) c -*y^* i 

v ; ^ ml 2 

m=0 

Here is a quick "proof" using the central limit theorem (this kind of 
argument goes back to Mark Kac [3]): Let X ly . . . ,X k be independent 
Poisson variables with parameter 1 (so having mean and variance 1). 
Let Sfc = X\ + • — \- Xk be their sum. By the Central Limit Theorem, 
s ^ k is asymptotically normal, in particular as k — > 00, 

Pro b(^^ < 0) ^ - 

that is 

Prob(,S fe < k) - i 

However since is Poisson with parameter fc, we have 

a- 

/ in 

Prob{S k <k)= e~ k — - 
^— ' m! 

m=0 

which gives (jI3J) . 

A. 2. The supremum of cusp forms. 

Proposition 3. Let f be a cusp form of weight k for SL 2 (1i). Then 
uniformly for z in compact subsets ofW, 

y k \f{z)\ 2 , 

~VJT <<k 

where (f, f) is the Peterson inner product. 
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Proof. Consider the integral of y k \f(z)\ 2 over the "Siegel set" 

6 := {z — x + iy : < x < 1, y > I/A-k} 

This integral converges since / is exponentially decreasing in the cusp. 
Since the Siegel set (5 is contained in a fixed number of translates of 
the standard fundamental domain T = {z = x + iy : \z\ > 1, \x\ < 1/2} 
we have 

\f(z)\V^ «</,/) 

s y 

On the other hand, using the Fourier expansion f(z) = ^ n>1 a(n)e 2nmz 
and Parseval, we find 



e 



\f(z)\ 2 y k ^ 

y2 



E 

n>l 



a[n) 



-A*ny y k-2 dy 



1/47T 



Thus we find that 



(14) 



E 



V^f r(*-i,n) 

n>l v ' 



^- 1 T(k-l,n)«(f,f) 



n> 



J (47m 



To estimate y k \f{z) | 2 , use the Fourier expansion and Cauchy-Schwartz: 



V*I/(*)| 2 = V* 



E«( 



n>l 



a(n)p 



< E 

n>l v 



jT(k-l,n)xy k J2^ 



■Airny 



(47m 



,fc-i 



n>l 



r(Jfe - l,n) 



For the first sum, use (|14|) to bound it by (/, /) . Thus we find 



y fc l/(*)F 
(fJ) 



«y fc E' 



k-l 



n>l 



T(k-l,n) 



To bound the sum, we split it up into "small" and "large" n's. For 
"small" n, that is those satisfying n < k, we can give an upper bound 
as follows: First use T(k — 1, n) > T(k — l,k), for n < k. This gives 



n<k 



-4nny 



(47TO 



.jfe-l 



< 



y 



T(k-l,n) ~ T(k-l,k) 

y 



hk) E' 



n<k 
oo 



< 



T(k-l,k) 



hk) E 



-4-7rnj/ 



,— iwny 



(4:7m) 



k-l 



(inny) 



k-l 



n=l 
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We have 

00 poo 

V e- A7Tny (Anny) k - 1 < / e~^ xy {ATixy) k - l dx + 2e~ {k - 1 \k - l) fe_1 

n=l ^° 

Thus we find 

^kST^ny yT{k) 1 1 

r(A;-l,n) T(A;-l,fc) l |/ ^ J 



n<A; 

which by ffTTj) is 



For the sum over the remaining n's, use (II U|) to get 

r(A;-l,n) > e" n n fc - 2 

which gives 

* V- -4*nv (47rn) fc - 1 < . 47tni/ (47rw) fc - 1 

n>k v ' y n>fc 



< (4vry) fc ^e- n(4 ^- 1 



n>k 



, A xk (k + l) e -(*+D(4^-i) 
« ^ (i- e -(^-i))2 



Anye k 



which, for fixed y > y/3/2, is negligible. □ 
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